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1  General  Information 

1.1  Status  of  Effort 

The  research  program  was  originally  planned  for  a  three-year  period.  However,  because  of 
the  new  appointment  of  the  principal  investigator  at  MIT,  the  program  was  divided  into  two 
stages.  The  first  stage,  carried  out  completely  at  the  Iowa  State  University,  produced  the 
following  major  results. 

•  A  study  of  dynamic  behaviour  of  rate  limiters  allowed  for  the  first  time  to  devel- 
ope  methods  of  absolute  stability  analysis  that  are  capable  of  discriminating  between 
systems  with  quasiconcave  nonlinearities  (such  as  the  saturation)  and  quasiconvex  non- 
linearities  (such  as  the  deadzone). 

•  A  software  package  for  rigorous  analysis  of  stability  and  performance  of  nonlinear /time- 
var)dng/uncertain  systems  was  developed,  based  on  the  concept  of  Integral  Quadratic 
Constraints. 

•  A  new  method  of  gain  scheduling  for  linear  feedback  design  under  the  condition  of 
control  saturation  has  been  developed,  with  guaranteed  stability  margins  rigorously 
proven. 

•  A  unique  technique  for  rigorous  analysis  of  global  and  semiglobal  stability  of  self¬ 
oscillations  in  higher  order  relay  feedback  systems  has  been  developed. 

1.2  Personnel  Supported 

Assistant  Professor  A.  Megretski,  Graduate  Students  A.  Cygankov  and  C.  Gunadi 

1.3  Publications  (submitted/accepted/published) 

1.  F.  Blanchini,  A.  Megi;etski,  Robust  State  Feedback  Control  of  LTV  Systems:  Nonlinear 
Is  Better  Than  Linear.  Submitted  to  IEEE  Transactions  on  Automatic  Control. 

2.  T.  Georgiou,  M.  Khammash,  and  A.  Megretski,  On  a  Large-Gain  Theorem.  Submitted 
to  Systems  and  Control  Letters. 

3.  A.  Megretski,  A.  Rantzer,  System  Analysis  via  Integral  Quadratic  Constraints.  Pub¬ 
lished  in  IEEE  Transactions  on  Automatic  Control. 
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4.  A.  Megretski,  Integral  Quadratic  Constraints  in  the  Design  of  Gain  Scheduled  Robust 
Controllers,  Published  in  the  Proceedings  of  CDC96. 

5.  A.  Megretski,  A.  Rantzer,  Integral  Quadratic  Constraints,  Part  I/Part  11.  Submitted 
to  Systems  and  Control  Letters. 

6.  A.  Megretski,  Global  stability  of  oscillations  induced  by  a  relay  feedback.  Proceedings 
of  the  1996  IFAC  World  Congress. 

7.  A.  Megretski,  L2  BIBO  output  feedback  stabilization  with  saturated  control.  Proceed¬ 
ings  of  the  1996  IFAC  World  Congress. 

1.4  Interactions/Transitions 

The  results  of  the  supported  research  have  been  communicated  in  presentations  at  the  1996 
Control  and  Decision  Conference,  the  1996  IFAC  World  Congress,  1996  MTNS  Conference, 
and  the  1997  European  Control  Conference.  The  developed  IQC  software  has  been  provided 
for  evaluation  at  the  Wright-Patterson  AFB,  and  at  Boeing  Inc.,  for  the  purpose  of  using  it 
in  the  analysis  of  gain  scheduled  flight  control  systems. 

2  Research  Accomplishments:  Technical  Details 

2.1  Integral  Quadratic  Constraints  based  on  lower  order  system  analysis. 

Integral  Quadratic  Constraints  (IQC)  are  used  to  describe  frequency-domain  input-output 
relations  in  common  nonlinear  and  time-varying  blocks.  Essentially,  in  the  IQC-based  analy¬ 
sis,  deterministic  relations  in  nonlinear  blocks  are  replaced  by  less  detailed  uncertain  models 
posessing  a  special  structure  which  allows  easy  analysis  via  convex  optimization.  The  quality 
of  such  analysis  depends  on  how  close  is  the  match  between  the  uncertain  IQC  model  and 
the  original  nonlinear  system  model.  While  in  the  past  a  number  of  IQC  models  has  ben 
derived  for  some  basic  classes  of  nonlinear  and  time-varying  blocks,  the  desired  quality  of 
approximation  has  not  been  achieved  yet  for  several  important  cases. 

In  particular,  the  classical  positivity /passivity  methods  did  not  provide  an  IQC  model  for 
the  saturation  nonlinearity,  which  would  distinguish  between  deadzone  and  saturation.  As  a 
result,  stability  analysis  of  systems  with  pure  integrators  controlled  via  a  saturated  feedback 
was  essentially  not  possible  using  the  classical  positivity/absolute  stability  framework. 
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Similarly,  the  analysis  of  gain  scheduled  control  algorithms  via  constant  D  —  G  scaling 
is  so  conservative  that  it  does  not  allow  to  account  for  any  positive  effects  of  nonlinear  gain 
scheduling. 

A  promissing  approach  to  derivation  of  advanced  IQC  models  is  described  below.  The 
approach  was  applied  to  obtain  less  conservative  stability  conditions  for  the  systems  with  pure 
integrators  controlled  via  saturated  feedback,  systems  with  backlash  hysteresis,  nonlinear 
gain  scheduled  systems,  relay/dry  friction  systems,  etc. 

Let  Aq  be  a  nonlinear  block  (in  many  cases,  Aq  is  a  memoryless  nonlinearity,  such  as 
relay,  controlled  gain,  saturation,  etc.,  or  a  nonlinearity  with  “discrete”  memory,  such  as  the 
hysteresis  nonlinearities).  We  consider  an  imaginary  experiment:  testing  Aq  in  a  feedback 
loop  with  a  low  order  (one  or  two  states)  linear  time  invariant  system  Go  (see  Figure  1). 
The  “results  of  the  experiment”  are  written  in  the  form  of  IQC  describing  the  input-output 


Figure  1:  Testing  Aq  in  a  low-order  feedback 

relation  between  /  and  z.  These  IQC  can  be  obtained  in  two  ways: 

•  re-writing  the  known  IQC  describing  the  relation  between  y  and  ^  as  IQC’s  describing 
the  relation  between  /  and  z; 

•  applying  advanced  system  analysis  tools  (non-quadratic  Lyapunov  functions,  optimal 
control  etc.)  in  the  analysis  of  the  closed-loop  system  in  Figure  1). 

Now,  when  the  problem  of  analysis  of  a  feedback  interconnection  of  Aq  with  other  nonlinear 
systems  is  considered,  the  interconnection  is  being  re-arranged  in  order  to  include  explicitly 
the  dynamic  block  on  Figure  1  (see  Figure  2).  There  are  two  sorts  of  benefits  from  the 
proposed  approach.  First,  new  important  IQC’s,  not  known  before,  can  be  obtained.  In  par¬ 
ticular,  an  exact  calculation  of  the  worst  case  gain  /  ->  z  in  the  feedback  interconnection 
of  a  saturation  and  a  pure  integrator 

y  =  sat(/  -y),  z  =  y  +  sat(/  -  y). 
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Figure  2:  Using  the  IQC  for  Gq/Aq  interconnection 

results  in  the  gain  value  7  =  \/2,  which  yields  perhaps  the  first  IQC  capable  of  distinguish¬ 
ing  between  the  saturation  nonlinearity  and  the  deadzone.  This  IQC  can  be  used  in  testing 
stability  of  systems  which  are  not  exponentially  stable  using  the  standard  gain  analy¬ 
sis.  Similar  results  are  being  developed  for  the  feedback  interconnection  of  a  double  pure 
integrator  with  a  saturated  feedback,  for  the  multiplication  nonlinearity  (important  in  the 
analysis  of  gain  scheduling),  etc. 

The  second  benefit  of  the  low-order  nonlinearity  testing  is  that  it  makes  possible  stability 
analysis  of  systems  with  bad-behaving  nonlinearities,  with  applications  to  the  nonlinearities 
such  as  the  ideal  relay,  dry  friction,  and  hysteresis. 

2.2  Special  IQC  for  semiconcave  nonlinearities 

The  main  technical  issue  in  this  subsection  is  validation  of  a  set  of  IQC  relating  signals  2 
and  X  in  the  system 

x(t)  =  f  (f){z{r))dr,  (1) 

Jo 

where  <f>  is  a,  “saturation-like”  memoryless  nonlinearity.  Though  system  (1),  because  of 
its  instability,  does  not  fit  well  within  the  standard  IQC  analysis  framework,  the  results 
originally  derived  for  (1)  can  be  easily  transformed  into  a  set  of  IQC  for  an  “encapsulated” 
rate  limiter  system,  defined  by 

x{t)  =  ^{v{t)  —  x{t)),  w{t)  =  x{t)  +  <f>{v{t)  —  x(t)),  x(0)  =  0  (2) 

(see  Figure  3).  For  example,  it  will  follow  from  the  main  result  that  the  gain  “from  v  to  w” 
in  system  (2)  does  not  exceed  y/2. 
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Note  that  while  the  gain  is  exactly  \/2  for  (^(y)  =  sat(y),  replacing  (f){y)  =  sat(y)  by  its 
linearization  at  zero,  ^(y)  =  y,  yields  an  identity  system  w  =  v  (the  induced  L2  gain  equals 
1),  while  replacing  ^(y)  =  sat(y)  by  <^(y)  =  dzn(y)  results  in  an  infinite  Lj  gain. 

X 


1/s 


s  +  1 


w 


Figure  3:  Encapsulation  of  a  rate  limiter 

The  IQC  result  has  broad  applications  in  the  analysis  of  more  complex  systems  (higher 
order,  other  nonlinearities,  time-variance,  uncertainty)  that  include  the  feedback  intercon¬ 
nection  (2)  as  a  subsystem.  Generally,  the  results  of  this  paper  can  be  applied  to  any  system 
of  the  form 


Vi 

V2 


=  Go 


AiVi 

<f>{v2) 


Go  = 


Gii  (s -h  1)((J12  —  l)/s 
(?2i  (s  -f  1)^22/ S 


(3) 


where  <j)  is  the  same  as  in  (2),  Ai  represents  other  nonlinearities/uncertainties  in  the  system, 
and  Gij  are  stable  proper  transfer  matrices.  While  system  (3)  is  not  given  in  the  standard 
IQC  analysis  format,  it  can  be  reduced,  via  a  simple  feedback  loop  transformation,  to  a 
standard  feedback  interconnection  of  a  stable  LTI  plant  with  a  structured  “uncertainty” 
which  consists  of  blocks  A  and  (2) 


Vi 

V2 


=  G 


Aivi 


G 


Gn  Gi2 
G21  G22 


(4) 


where  A^  is  the  operator  v  w  defined  by  (2)  (see  Figure  4). 


Figure  4;  Loop  transformation  for  encapsulation 


The  results  of  this  section  hold  for  a  large  class  of  semiconcave  functions  </>.  The  definition 
of  a  semiconcave  function  summarizes  those  features  of  the  ideal  saturation  nonlinearity 
which  are  essential  in  proving  the  \/2-gain  result  and  its  generalizations. 
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Definition  A  monotonically  non-decreasing  odd  function  ^  :  R  — >  R  is  called  semiconcave 
if  <l){z)  =  <j>{z)/z  is  monotonically  non-increasing  over  the  interval  z  G  (0,  oo).The  set  of  all 
semiconcave  functions  will  be  denoted  by  SC. 

It  is  easy  to  see  that  a  semiconvcave  function  (f)  =  (f){z)  is  differentiable  at  z  =  0  iff  ^{z)lz 
is  bounded  on  the  interval  (0,  oo),  in  which  case 

(^(0)  =  sup  (5) 

z>0  Z 

For  convenience,  we  will  consider  (5)  as  a  definition  of  (^(0)  in  the  case  when  the  right  side 
in  (5)  is  infinity.  Note  also  that  0(0)  =  0  would  imply  0  =  0,  in  which  case  system  (1)  is 
trivial. 

Theorem  2.1  Let  <{>  be  a  semiconcave  function  with  0(0)  =  K ,  0  <  K  <  oo,  b  ^  "R.  The 
inequality 

{2\z  +  bx\^-\z\'^-z(f>{z)/K}dt>0  (6) 

holds  for  any  x,z  G  L2  satisfying  relation  (1).  Moreover,  if  b>0  then 

^  {2|z  +  top  -  Izp  -  z<t){z)/K}  dt>0  (7) 

for  allT  >0. 

Inequality  (6)  can  be  considered  as  a  family  of  IQC  describing  the  relation  between  z 
and  0(z).  Since  the  corresponding  ai  is  defined  by  an  unstable  system  L  (contains  a  pure 
integrator),  it  is  difficult  to  use  (6)  directly  in  the  analysis  of  systems  involving  semiconcave 
nonlinearities.  To  resolve  this  problem,  (6)  is  re-written  as  a  set  of  IQC  describing  the 
“encapsulated”  rate  limiter  (2)  (see  Figure  3). 

Theorem  2.2  Let  f)  be  a  semiconcave  function  with  0(0)  =  K,  0  <  K  <  00,  6  G  R,  a  >  0. 
Define  as  a  system  v  w,  where 

x  =  (f>{v  —  x),  w  =  ax  +  <f>{v  —  x),  x(0)  =  0.  (8) 

The  system  is  stable,  and  the  inequality 

2\\v  +  bx\\‘^  >  llu  —  a;|p  +  {v  —  x,w  —  ax)/K  >  0  (9) 

holds  for  any  v  G  L^. 
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Proof  of  Theorem  2.2.  For -6  =  1,  z  =  v  —  x,  inequality  (7)  implies  ||Pt(w  —  a:)||  <  2||Ptu1|. 
Hence  system  AJ  is  stable.  In  particular,  x,w  E  L2  whenever  v  E  L2.  Now  (6)  implies  (9) 
for  any  6  G  R.  ■ 

Note  that  x  —  GaW  where  Go(s)  =  l/(s  +  a),  i.e.  (9)  are  true  IQC  relating  v  and  w.  If  ^ 
is  not  a  semiconcave  function,  these  IQC  are  generally  not  valid.  For  example.  Theorem  2.2 
with  K  —  l,a  =  l,b  =  —0.5  implies 

||t;|p  >  0.5||x||^  +  {v  —  X,  (f){v  —  x))  >  0.5||x  +  <l){v  —  x)|p, 

where  it  was  used  in  the  derivation  that 

(x,  (j){v  —  x))  =  {x,  x)  =  0, 

and 

{v  —  X,  <f>{v  —  x))  >  \\(j){v  -  x)||^. 

Therefore,  ||A^||  <  \/2  when  ^  is  semiconcave  and  0(0)  <  1.  On  the  other  hand,  it  is  easy 
to  see  that  replacing  0  by  dzn  yields  an  unstable  system  A^  (infinite  L2  gain). 

There  are  other  IQC  known  to  describe  the  relation  between  z  and  (j){z).  The  criterion 
by  Zames  and  Falb  (SIAM  J.  of  Control,  6(1):89-108, 1968)  states  that 

{Kz  -  0(^),  Hiz)  +  {H2{Kz  -  <t>{z)\  z)  >  0,  (10) 

where  Hi,H2  are  LTI  systems  such  that 

poo 

Hiis)  =  Di  +  Ci{sl  -  Di  >  /  \Cie^'%\ dt,  (11) 

Jo 

and  0  G  [0,  AT].  These  IQC  can  be  re-written  as  input/output  relations  for  w  =  AJw: 

{K{v  -  x)  -  {v  -  aw),Hi{v  -x))  +  {H2{K{v  -x)  -  (v  -  aw)),v  -  x)  >  0,  (12) 

where,  as  before,  x  =  GaW. 
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